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We present a detailed analysis of bi-partite entanglement in the non-Abelian Moore-Read fractional quantum 
Hall state of bosons and fermions on the torus. In particular, we show that the entanglement spectra can be 
decomposed into intricate combinations of different sectors of the conformal field theory describing the edge 
physics, and that the edge level counting and tower structure can be microscopically understood by considering 
the vicinity of the thin-torus limit. We also find that the boundary entropy density of the Moore-Read state is 
markedly higher than in the Laughlin states investigated so far. Despite the torus geometry being somewhat more 
involved than in the sphere geometry, our analysis and insights may prove useful when adopting entanglement 
probes to other systems that are more easily studied with periodic boundary conditions, such as fractional Chern 
insulators and lattice problems in general. 



PACS numbers: 73.43.Cd, 71.10.Pm, 03.67.-a 

I. INTRODUCTION 

Quantum correlations give rise to many exotic phases of 
matter that cannot be characterized in terms of traditional 
concepts, such as local order parameters and symmetry. Re- 
cently, tools from the field of quantum information (QI) have 
been used to quantify such correlations'. Of special inter- 
est among the applications are systems in which more tradi- 
tional condensed-matter methods are of limited use, for exam- 
ple topologically ordered matter*^. Fractional quantum Hall 
(FQH) states stand out as experimentally verified topologi- 
cally ordered phases driven by interactions, and their possi- 
ble applications in the context of quantum computation are 
of great current interest^. The microscopic understanding 
of these phases is mainly based on ad hoc, albeit brilliant, 
guess work''"^^ and numerical wave-function overlap calcula- 
tions in small systems. A fundamental problem with using 
wave-function overlaps as a probe is, however, that it nec- 
essarily vanishes in the thermodynamic limit (for any real- 
istic interaction). Recently, it has been realized that (bi- 
partite) entanglement measures, most saliently the von Neu- 
mann entropy'" " and the entanglement spectrum'- can pro- 
vide valuable insights into these states — in principle even in 
the thermodynamic limit. 

In this work, we focus our attention on entanglement in the 
archetypical non-Abelian FQH state, namely the Moore-Read 
state^, which has received a tremendous amount of attention 
recently as a potential platform for topological quantum com- 
putation. Previous theoretical studies'^ have accumulated ev- 
idence that the ground state of the two-dimensional electron 
gas at the Landau level filling fraction ly — 5/2 is well de- 
scribed by the Moore-Read state, which may be thought of 
as paired composite fermions and has quasiparticles possess- 
ing fractional charge ±e/4 and obeying non-Abelian braid 
statistics"". In recent experiments, both the fractional charge 




FIG. 1. (Color online) The torus setup (a) compared with the orbital 
partition on the sphere (b). The dark lines indicate the centers of 
the single-particle states and the differently shaded regions denote 
the approximate spatial partitioning corresponding to the half-block 
orbital partitioning. Red arrows represent the artificial edge states 
induced by splitting the system into A and B. 

and non-Abelian braid statistics have been claimed''*, but 
the interpretations of the experiments is still under debate'"". 
Another possible host of the Moore-Read state is the Bose- 
Einstein condensate under rapid rotation, in which the bosonic 
state at = 1 is particularly promising"'. However, the ex- 
perimental realization of the bosonic FQHE is extremely chal- 
lenging, although some strategies to overcome the difficulties 
have been proposed'^. 

To study bi-partite entanglement, we (artificially) divide a 
system into two parts A and B (Fig. 1). In a tensor product 
Hilbert space, H ~ Ha 'E)'Hb, any pure state can be 

decomposed using the Schmidt decomposition'**, 

i 

where the states lipf) (lipf)) form an orthonormal basis for 
the subsystem A (B) and the entanglement "energies" > 
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are related to the eigenvalues, Xi, of the reduced density 
matrix, pA = tiB \ '^)ab ab(*|, of A as Xi = e"?'. 

For topologically ordered states in two dimensions, the en- 
tanglement entropy contains topological information about 
the state: Sa = -tr[p^lnpA] = -J^iX^^^^i = 
J2i ^jC^^S is expected to scale as 

SA~aL-nj + 0{l/L), 

where L is the (total) block boundary length, n is the number 
of disconnected boundaries, and 7 characterizes the topologi- 
cal field theory describing the state'" " . 

Li and Haldane'- realized that the full so-called entangle- 
ment spectrum (ES), {£,i}, contains much more information 
than entanglement entropy. In particular, when plotted against 
the natural quantum numbers of the system, it shows a remark- 
able similarity with the conformal field theory (CFT) describ- 
ing the chiral edge states of the FQH states. 

To make practical use of the entanglement concepts, it is 
instrumental to find a protocol with which the theoretical 
ideas can be (numerically) tested in realistic circumstances. 
The most widely used concept of partitioning the system in 
terms of the single-particle orbitals was introduced by Haque, 
Zozulya, and Schoutens"" in their study of the topological en- 
tanglement entropy in Laughlin states on the sphere. A numer- 
ical determination of 7 (and a) in realistic circumstances re- 
quires information about Sa for a number of different bound- 
ary lengths, L. Because of its technical simplicity, early at- 
tempts to obtain the entropy scaling in FQH states focused 
on the sphere geometry^". However, as recently demonstrated 
for Abelian FQH states, a substantially better finite size scal- 
ing can be obtained on the torus where the boundary length 
can be varied continuously by varying the aspect ratio^' [cf. 
Fig. 1(a) and I (b)]. (The idea of obtaining entanglement 
entropy scaling through varying discrete torus circumferences 
was also used in Ref. 22 for the dimer model on the triangular 
lattice.) Importantly, this extra degree of freedom available on 
the torus also provides a handle on when the extrapolations 
needed to extract 7 can be trusted (and when they cannot). 

With a few very recent exceptions^^*"^^, the efforts made in 
the study of the ES in FQH states are also numerical 
In addition to these works, there has been a large number of 
recent studies extending the range of applicability of the ES 
to an increasing number of physical systems'***. The studies 
of the ES in FQH states have focused predominantly on the 
sphere geometry. In this case, there is a genuine benefit with 
this choice since it amounts to probing the physics of a single 
FQH edge while the natural partition on the torus corresponds 
to two oppositely oriented edges [cf. the red arrows in Fig. 
1(a) and 1 (b), respectively]. A benefit with the torus setup 
is, however, that one can continuously connect to the exactly 
solvable thin-torus limit^** from which many of the properties 
of the ES can be understood microscopically^*. 

On the sphere one finds that the ES has a chiral structure'^ 
that is intimately related to the squeezing rule of model FQH 
states'*" that holds on genus-0 manifolds. The structure of the 
squeezed configurations also provides physical insight similar 
to what is possible in the thin-torus limit. While the squeezing 
rule does not hold on the torus (genus- 1), the ES can never- 



theless be described by combining two edge spectra, as was 
shown in Ref. 28 for the Laughlin state. 

In spite of the technical difficulties involving two separate 
edges, these issues are worth dealing with, in particular since 
there are many physical systems of great interest that are only 
approachable using periodic boundary conditions. Specifi- 
cally, regular two-dimensional lattices do not admit generally 
a defect-free embedding onto the sphere (because of their dif- 
ferent Euler characteristics). In particular, the recently pro- 
posed fractional Chern insulators'*''"*^ appear to belong to this 
category. 

The two-edge picture on the torus is reportedly'''' '*^ difficult 
to extend to non-Abelian FQH states due to their non-trivial 
ground-state degeneracies, which do not result from simple 
center-of-mass translations as in the Abelian case. Thus, it is 
not a priori clear how to choose the ground state, in 
(1) (or alternatively, how to define the density matrix of the 
full system AU B) out of this degenerate set. Note that the 
issue of degenerate ground states does not occur in the sphere 
case in which the model states are unique maximal density 
zero modes of their respective parent Hamiltonians. 

Here, we adopt a very simple and natural choice for the set 
of \'i')AB and show that a similar, but significantly richer, two- 
edge picture also holds true for the ES of non-Abelian FQH 
states on the torus. Specifically, we disentangle the physics of 
the edge modes appearing in the entanglement spectra in each 
of the topologically distinct sectors of the Moore-Read state 
of both fermions and bosons. We find that, even for a given 
cut in one of the ground states, the resulting towers are gener- 
ated from combinations of different sectors of the underlying 
conformal field theory. 

We also carefully analyze the scaling of the von Neumann 
entropy in the various sectors of the Moore-Read state. We 
find that the total entropy as well as the area-law entropy 
density, a, can be estimated (in particular quite accurately in 
the case of bosons) while the extrapolation is too sensitive to 
faithfully determine the topological part, 7. 

The remainder of this article is organized as follows. In 
Section II, we introduce the physical model and the method 
we use to obtain the ground states and calculate the ES. In 
Section III, we analyze the ES from two distinct perspectives. 
On the one hand, we explain the ES as the combination of 
edge modes and discuss the quantitative relation in this com- 
bination. On the other hand, we use the thin-torus limit and 
a perturbation theory to illuminate the microscopic origin of 
the observed ES, including the counting rules in different edge 
sectors. Finally, we discuss the entanglement entropy in Sec- 
tion IV. 



II. MODEL AND METHOD 

We study a two-dimensional A^-boson (fermion) system 
subject to a perpendicular magnetic field on a torus with peri- 
ods Li and L2 in the x and y directions. The full symmetry 
analysis of this system was first provided by Haldane'*^ — here 
we use a convenient representation thereof. Periodic boundary 
conditions require that L1L2 — 2ttNs (in units of the mag- 



3 



netic length) where Ng is the (integer) number of magnetic 
flux quanta (the number of vortices for rotating Bose-Einstein 
condensates). We choose a basis of normahzed single-particle 
lowest Landau level (LLL) wave functions as 
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lii?^+nL^)x-(y+nL^+?^f/2] 



(2) 



where j = 0, 1, 2, A^^ — 1 can be understood as the 
single-particle momentum in units of 2tt/Li. Because ijjj is 
centered along the line y — — 27rj/Li, the whole system can 
be divided into Ng orbitals that are spatially localized in the 

direction (but delocalized in the .t— direction). There are 
two translation operators, T^, a ~ 1.2, that commute with 
the Hamiltonian H (and any translational invariant operator); 
they obey T1T2 = e^'^'^/^-TaTi, and operators have eigen- 
values e^'^'-^'^/^^ifa = 0,...,Ns - 1. Ti corresponds to 
x-translations and Ki — J^iLiJi (mod Ng) is the total x- 
momentum in units of 2tt/Li. T2 translates a many -body state 
one lattice constant L2/NS — 2tt/Li in the y-direction and 
increases Ki by N . At filling factor v ~ p/q (with p and q 
co-prime), r| commutes with Ti, and T2 (fc = 0, 1, . . . g — 1) 
generate q degenerate orthogonal states, which have different 
Ki. This is the g-fold center of mass degeneracy common to 
all eigenstates of a translational invariant operator in a Lan- 
dau level. Thus, the energy eigenstates are naturally labeled 
by a two-dimensional vector — 0, Ng/q — 1, where 



e'^^'^'^^i/^s is the T|-eigenvalue. 

We use exact diagonalization to obtain the Moore-Read 
states, which are zero-energy ground states of certain three- 
body Hamiltonians (see Appendix A), in the orbital basis. The 
Moore-Read states are non-Abelian states, for which the de- 
generacy on the torus is enhanced (in this case by a factor 3) 
compared to the g-fold degeneracy discussed above. It is read- 
ily seen from the thin-torus configurations (the ground states 
as Li 0) that they are not simply the translations of each 
other"*"* (see below). To extract the ES, we choose the ground 
states as eigenstates of Ti and r| and bipartition the system 
into blocks A and B, which consist of I a consecutive orbitals 
and the remaining Ns — Ia orbitals, respectively. We label 



every ES level by the particle number Na = SjgA "^j 
the total momentum Ka — SjgA (™od N^) in block A, 
where Uj is the particle number on the orbital j. (In this work, 
we present data only for the case in which I a = Ns/2.) 

To understand the ES, it is essential to understand what the 
partitioning of the state looks like in the thin-torus limit. For 
the bosonic case, there are three different thin-torus patterns 
leading to the following partitions (for N = Ng = 16): 



1111|11111111|1111 
020210202020210202 ± 202012020202012020. 



(3) 



For the fermionic case, there are six different thin torus pat- 
terns and the following partitions (for N = 16, Ng — 32): 



01010101|0101010101010101|01010101 

10101010|1010101010101010|10101010 
01100110|0110011001100110|01100110± 10011001|1001100110011001|10011001 
1100110011100110011001100111001100 ±0011001110011001100110011100110011. 



(4) 



The bold block is our subsystem A. For bosons in (3), 
we have two qualitatively different cuts: 11|11---11|11 and 
02|02 • • • 02|02 (20|20 • • • 20|20 gives a mirror image of this). 
For fermions in (4), we have four qualitatively different cuts: 
01|01---01|01 (10|10---10|10), 0110I0110-- -011010110, 
1001|1001---1001|1001, and 1100|1100 • • • 1100|1100 
(0011|0011---0011|0011). 

We stress that, as long as the edges are sufficiently well 
separated, one can understand the entanglement in terms of 
two non-interacting edges whose details depend on the local 
environment around the cuts^*^. This holds true also for the 
states that are connected to a thin-torus configuration which is 
a linear superposition of two individual terms — in these cases 
the ES is composed of two shifted and superimposed mirror 
images corresponding to the ES of a single term respectively. 

Our procedure is different from that in Refs. 33 and 42 
where the authors calculate the ES via a mixed state den- 
sity matrix of the form p = \ X^iLi \^^)ab ab{^^\ where 
denote d degenerate ground states. With this recipe 
one finds that the ES corresponds to the superimposed ES of 



all the d thin-torus patterns. For the entanglement entropy, 
such a mixed state prescription essentially shifts Sa{L) by 
a constant and would thus result in a shifted prediction for 
the topological contribution, 7. In the case of Abelian states, 
it turns out that averaging the entropies (rather than the den- 
sity matrices) over the different sectors, or equivalently over 
the possible translations of the region A, significantly reduces 
finite-size corrections and yields results in excellent agree- 
ment with theory-'. We note that the mixed-state prescription 
shifts the entropies of Abelian states by a constant value. In d, 
and would thus lead to a topological entropy different from the 
theoretical predictions for the spatial (as opposed to orbital) 
cut — in fact, it would lead to 7 = 0. For non-Abelian states, 
it is not yet settled which orbital basis prescription would lead 
to the same topological entropy as for the spatial cut. 
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FIG. 2. (color online) The edge modes of the environment in the left 
panel (black dots) and the environment in the middle panel (green 
dots) can combine to form a tower in the right panel (blue crosses). 
The relation in Eq. (5) is shown by the parallelogram. The edge 
mode at Afc = — 1 pointed by the solid black arrow and the edge 
mode at Afc — 1 pointed by the solid green arrow can generate the 
level at Afc = pointed by the solid blue arrow. This data comes 
from the ES of bosons in the 1 1 sector (see Fig. 5). 



III. ENTANGLEMENT SPECTRA: TWO-EDGE PICTURE 
AND THIN-TORUS ANALYSIS 

The most prominent Na sectors of the ES of the Moore- 
Read state for iV 16 are displayed in Fig. 5 (v — 1 bosons) 
and Fig. 6 (i/ = 1/2 fermions). The gross features of the 
Moore-Read ES on the torus are very similar to that of the 
Laughlin state — in both cases, multiple towers are formed"**. 
In this section, we analyze the ES from two different perspec- 



tives: We explain the tower structure in terms of combinations 
of edge modes and highlight intriguing relations between the 
ES levels within the towers as well as between the levels in 
different particle number sectors. Moreover, we use the ex- 
actly solvable thin-torus limit and perturbation theory to un- 
derstand the formation of various edge environments and tow- 
ers. 

The observed towers in the numerical ES can be reproduced 
by first assigning the edge modes of individual edge environ- 
ments and then combining them appropriately. The number 
of independent edge modes at momentum Afc in an edge en- 
vironment is determined by the underlying edge theory. The 
edge theory of the Moore-Read state is richer than that of the 
Laughlin state and contains a free boson branch as well as a 
Majorana fermion branch"**". The details are recapitulated in 
Appendix B for completeness. It is important to note that 
there are different sectors of the edge theory and that they 
come with different predictions for the counting of states as 
a function of momentum. This is reflected in our numerically 
obtained ES, where we observe the edge environments with 
different counting rules. It is interesting to see that two edge 
environments with different counting rules can also combine 
to form a tower 

There are intriguing quantitative relations in the combina- 
tion of edge modes as first pointed out for the Laughlin state 
in Ref. 28. An explicit example of how two edges, with dif- 
ferent dispersion, add up to a tower is given in Fig. 2. More 
generally, each edge mode can be labeled by three parameters: 
the edge environment X to which it belongs, its momentum 
shift Aki compared with the bottom mode of the environment 
X, and the change of the subsystem particle number AN^ in 
the environment X compared with the thin-torus state. Two 
edge modes with entanglement energy £,{X, Aki, ^^a) ^'^'^ 
^(3^, Akj, AN^), respectively (here we assume Afc^ < and 
Afcj > 0), combine to form a level in the Xy tower with 
entanglement energy 



X 



Akj, AN A 



ANj) 



Ah, AN^) + ^k^ANX) ~ -[ax, 0, AN^) + ay, O, ATV^)]. (5) 



r 



The validity of the two-edge picture is insensitive to the cir- 
cumference Li as long as the edges are sufficiently well sep- 
arated from each other, i.e. given that d ^ L2/2 ~ ttNs/Li 
is large enough, which is equivalent to small enough Li for 
a given system size. This is illustrated in Fig. 3, where the 
breakdown of the two-edge picture is signaled for the larger 
Li values, which is indeed a confirmation of the fact that the 
decomposition of the entire ES into a combination of edge 
modes is a highly non-trivial fact. Note that this breakdown 
occurs despite the fact that the numerically exact Moore-Read 
state is obtained for all Li. 

The relative pseudo-energies of the assigned single-edge 
modes depend smoothly on the torus thickness Li, to some 
degree even after the two-edge prediction breaks down. For 



a given Li, the edge levels correspond well to the single- 
edge levels extracted from the ES on a sphere with a cor- 
responding length of the equator, as shown in Fig. 4. For 
large boundary lengths, the dispersion of a single edge be- 
comes non-monotonic — at least in the case of fermions [cf 
Fig. 4(b)], as can be inferred from the original data obtained 
by Li and Haldane'^. This does not imply that the two-edge 
picture will eventually break down, but does imply that the 
edge assignment becomes much more cumbersome at large 
Li as the Afc = levels no longer play the role of vacuum 
levels of each tower. Also, for this reason it is very useful to 
follow the evolution of the edge levels down to small Li where 
the dispersion is monotonic in order to eventually understand 
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FIG. 3. (color online) A plot of the main tower(s) of the ES in the 
0101 fermionic Moore-Read state for various Li (Na = 8,N = 
Ns/2 = 16). For small enough Li (in this case Li < 7 or so) the 
edges are well enough separated (d ~ L2/2 — nNs/Li) and the 
two-edge prediction (black crosses) reproduces the numerically ob- 
tained ES levels (blue squares). For larger Li the edges are spatially 
closer and the two-edge prediction gradually breaks down. 



the ES also at large Li . 

The adiabatic connection to the thin-torus (Li — > 0) limit 
also enables us to understand more detailed features of the 
ES by perturbing away from this solvable limit^**. The per- 
turbation theory is however hard to perform in a rigorous 
way, as was recently performed for the ES of one-dimensional 
models'*''. The reason for this is that the exponential behavior 
of the matrix elements implies that higher-order contributions 
from local terms come with amplitudes of the same order as 
longer-range terms contribute at lower orders. Nevertheless, 
many insights can be gained from a perturbative perspective, 
as we discuss below. 

It is instructive to divide the perturbations, which are three- 
particle hopping processes, into three different classes. In the 
first class, three particles belong to the same subsystem and 
none of them move across the edge. These processes do not 
qualitatively alter the entanglement between two subsystems. 
In the second class, two particles belong to one subsystem and 
one particle belongs to the other, but still none of them moves 
across the edge. In the third class, some of the particles move 
across the edge. As we show below, the processes in the sec- 
ond class are responsible for generating new levels within a 
tower, and those in the third class lead to levels in new tow- 
ers stemming from new edge environments. The entire ES of 
the Moore-Read state is built from successive combinations of 
many of the processes in each of these three classes. 

With the knowledge of the microscopic environment near a 
cut in the thin-torus limit, the counting of each edge follows 
from the exclusion principle that no more than two bosons 
(fermions) occupy two (four) adjacent orbitals. Similar ex- 
clusion rules are, in addition to the thin-torus limif*"*, also 



FIG. 4. (color online) The single-edge modes identified from the ES 
as a function of Li. (a) The edge modes corresponding to the 0|2 cut 
in the bosonic Moore-Read state for N — Ns = 12. (b) The edge 
modes corresponding to the OjO cut in the fermionic Moore-Read 
state for — 12, Ns = 24. The rectangles contain the single- 
edge ES levels in the spherical geometry here shifted: {^i} — > 
{^i + constant}, for the best comparison with the torus results at 
Li = 14 for bosons and Li = 16 for fermions. The number of flux 
quanta on the sphere N^'' is chosen as the integer nearest to / (27r^ ) 
so that the length of the equator of the sphere is nearly the same with 
Li. Here N = 12,A^|p = 10 for bosons and A^ = 8, Af = 13 
for fermions"***. One can see that in (b) the red dot is slightly lower 
than the black circle at Li — 13, meaning that a non-monotonic 
dispersion of the edge appears. 



showing up in related approaches'" ''^ such as the squeezing 
rules related to Jack polynomials and the patterns of zeros ap- 
proach. 



A. Bosons 

We now give a more detailed account of the Moore-Read 
ES in the case of bosons. 

We first consider the 1 1 sector and systematically explain 
the ES in this sector, which is shown in Fig. 5. The lowest ES 
level is found in the Ny^ — 8 sector at AKy^ = 0, correspond- 
ing to the thin-torus configuration 

1111|11111111|1111. 

At Li = 0, this is the only entanglement level. We call the 
edge environment 1111|1111 Ab, where the subscript, B, in- 
dicates the environment is for bosons. Here the subsystem on 
the left (right) side of this edge environment is A when we 
consider the right (left) edge of A. By definition AA'^" — 0. 
(In the following, when we discuss AN^ of an edge environ- 
ment X, we suppose the subsystem A is on the left side of 
X. If A is on the right side, one only needs to put a minus 
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FIG. 5. (color online) The ES of bosonic Moore-Read states in the 1 1 sector at Li — 5.5 (upper panels) and the 02+20 sector at ii = 6 (lower 
panels) for N = Ns — 16. The origin of AKa is chosen to match the Tao-Thouless state. The blue squares represent numerically obtained 
data. The assigned edge modes are labeled by dots with different colors corresponding to different edge environments we describe in the text. 
The combination of two identical edge environments is marked by the crosses with the color of that edge environment. The combination of two 
different edge environments is marked by the squares filled with two colors, the left (right) one of which corresponds to the edge environment 
on the left (right) edge of the subsystem A. Here we do not differentiate the edge environments V'^ to Xg from the edge environments Db to 
Ib because the former ones are just the mirror symmetries of the latter ones and have a momentum shift ±Ns/2. 



sign before the value that we give.) All other levels in the 
.4b .4b tower (y^B-4B denote the edge environments on the left 
and right edge of subsystem A) are generated from this level 
by the momentum-conserving hopping processes, which con- 
serve Nji. For example, a hopping process at the right edge 
of subsystem A, 1111|1111 1102|2011, gives the lowest 
level at AKa = 1. 

Some processes do not conserve Na', e.g., 1111|1111 ^■ 
1103|0111 or 1111|1111 1110|3011. We call the new 
kind of edge environment in this example Bb- It is clear that 
AiV^" — ±1. Two Bb edges create the B^Bb tower in the 
Na = 8 sector, whose dominant thin-torus configurations are: 

1103|01111103|0111 
1110|30111110|3011. 

Similarly, we can find another new edge environment. 
When applying a hopping process to the edge environment 
Bb: 1103|0111 ^ 1005|0011 or 1110|3011 ^ 1100|5001, 
we obtain the edge environment Cb with AA^^" = ±2. Two 
Cb edges create the CbCb tower in the A^^ — 8 sector, whose 
thin-torus configurations are: 

1005|00111005|0011 
1100|50011100|5001. 



The different edges can combine with each other to form 
towers in other Na sectors. For example, in the Na = 7 
sector we predict and observe AbBb, BbAb, BbCb and Cb.Bb 
towers. In the Na = 6 sector, we can observe another BbBb 
tower and ^Ib^b and CbAb towers. 

The .4b, Bb and Cb edges are sufficient to accurately re- 
produce the ES of the 11 sector up to ^ = 30 for Li — 5.5, 
as shown in the upper panels of Fig. 5. For larger Li more 
towers appear and can be explained along the same lines. 

Now we turn to the (asymmetric cut in the) 02H-20 sec- 
tor, for which the ES possesses more complicated structures 
than that in the 11 sector, as shown in Fig. 5. For simplic- 
ity, we start our analysis from only one term in the superpo- 
sition of the thin-torus configuration, for example, from the 
term 0202 102020202| 0202. The entire ES in the 02H-20 sec- 
tor is recovered by superposing two mirror images of the ES 
stemming from the single term. 

In the thin-torus limit, the only entanglement level at 
AKa = corresponds to the configuration 

0202|02020202|0202. 

We refer to this edge environment 0202 1 0202 as X>b- All 
other levels in the Vb'Db tower are generated from this level 



7 



by the momentum-conserving hopping processes, which con- 
serve Na- 

Through applying leading hopping processes, we can gen- 
erate all edge environments we have observed in the ES (at 

Li = 6): 





0202 10202, A A^J" = 0; 


£b 


0201 1 2102, ATV^" = -1 


Tb 


02001 4002, AA^J"* = -2 


Qb 


0104|0102, A7V^« = 1; 


T~Lb 


0006 10002, A A^;^» = 2; 


Xb 


0100|6001, AA^^" = -3 



These edges combine with each other to form towers in 
each Na sector For example, in the Na ~ 8 sector we predict 
and observe VbVb, Sb^b, ^b-^b, GbQb and HbHe towers, li 
the Na — 7 sector we find VbSb, SbJ^b, GbT^b, T-LbGb and 
J-bIb towers, and in the Na ~ 6 sector we find VbJ-b, Gb£b, 
Hb'Db and Eb^b towers. 

Similarly, we can start the analysis from the other term 
2020|20202020|2020 in the thin-torus configuration. We also 
predict and observe the six edge environments Vg to Zg, 
which are just the mirror symmetries of the edge environments 
Vb to Ib for the term 0202 1 02020202] 0202. For example, the 
edge environment I?g is identified as 2020 1 2020, which is the 
mirror symmetry of T>g. Moreover, the combination of edge 
environments to Ig can form towers in each Na sector. 
For example, in the Na = 7 sector we can observe the EgV^ 
tower as the mirror symmetry of the VbSb tower. 

We are also able to understand the counting rule of each 
edge environment from a simple exclusion rule in their thin- 
torus configuration. Here we take the edge environment Ab in 
the 1 1 sector and 2?b in the 02H-20 sector as examples. When 
analyzing the counting rule, we imagine the subsystem on the 
left (right) side of the edge environment as a quantum Hall 
system with an open right (left) edge, and then we move parti- 
cles to the orbitals with higher (lower) momentum to increase 
(decrease) the momentum of the system. Meanwhile, the gen- 
eralized exclusion rule'*" '*'* of the Moore-Read state, namely 
no more than two bosons (fermions) on two (four) consecu- 
tive orbits, should not be violated. Through the analysis, we 
can find that the counting rule of the edge environment Ab in 
the 1 1 sector is consistent with that of free bosons plus peri- 
odic Majorana fermions, while the counting rule of the edge 
environment Vb in the 02H-20 sector is consistent with that of 
free bosons plus antiperiodic Majorana fermions with an even 
F (see Appendix C). For some edge environments, only one 
side of it satisfies the generalized exclusion rule, for example 
the edge environment £b = 0201|2102 in the 02H-20 sector 
In this case, we only need to analyze the subsystem on its left 
side. Through analysis, we can predict the counting rules of 
some edge environments and compare them with the counting 
rules that we observe in our numerical data. In the 1 1 sector. 



we have 

^b: 1,2,4,8,14,--- ; [1,2,4,3]; 
Bb: 1,2,4,8,14,... ; [1,2,3,1]; 
Cb: 1,2,4, 8,14,... ; [1,2,1] 

and in the 02H-20 sector, we have 



■Db 


1,1,3,5,10,--- 


[1,1,3,3, 


£b 


1,2,4,7,13,--- 


[1,2,3,2, 


Tb 


1,1,3,5,10,--- 


[1,1,2,1]; 


Gb 


1,2,4,7,13,--- 


[1,2,2]; 


T~Lb 


1,1,3,5,10,--- 


[1,1]; 


Tb 


1,2,4,7,13,--- 


[1], 



where we, for each edge environment, first give the expected 
counting rule and then list the observed result from our nu- 
merical data, as shown in Fig. 5, in the brackets. For ex- 
ample, Ab ■ 1, 2, 4, 8, 14, • • • ; [1, 2, 4, 3] means that the ex- 
pected number of edge modes for edge environment Ab is 
1, 2, 4, 8, 14, • • • at Ak = 0, 1, 2, 3, 4, • • • , while the observed 
number of edge modes is 1, 2, 4, 3 at Ak = 0, 1, 2, 3. We see 
that the numerically observed count never exceeds the theo- 
retical expectations (that are derived for an infinite system). 
There are two reasons for this. First, our data are only numer- 
ically accurate up to some finite ^, and thus we count only the 
levels that are free of numerical noise. Second, the counting is 
truncated by the finite system size similar to the situation on 
the sphere'^"'* and should be expected in any geometry. 



B. Fermions 

The thin-torus and edge analysis of the fermion ES (Fig. 6) 
is entirely analogous to the boson case, and thus we only pro- 
vide a condensed exposition of the analysis here. Note how- 
ever, that the edge assignment would have been much trickier 
in the fermion case if we would have started out by consid- 
ering the large Li regime where the edge dispersion is non- 
monotonic'^. 

In the 0101 sector, we can observe three edge environments: 

Af: 01010101)01010101, AiV;^" = 0; 
Bp: 01010100)11100101, AiV^f = -1; 
Cp: 00011111)00010101, AiV^" = 1. 

Their counting rules are 

A: 1,2,4,8,14,... ; [1,2,4,4,1]; 
6f: 1,2,4, 8,14,... ;[1,2, 3,1]; 
Cf: 1,2,4,8,14,... ;[1,1], 

which are the same as those for the 111 bosonic state. The 
combinations of edge environments can form towers in each 
Na sector: y^p-^p, BpBp and CpCp towers in the Na = 8 
sector, ApBp and Cp^p towers in the Na — 7 sector, and 
CpBp tower in the Na = 6 sector. 
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FIG. 6. (color online) The ES of fermionic Moore-Read states in the 0101 sector at Li = 7 (upper panel), the 01 lO-l-lOOl sector at Li = 8 
(middle panel) and the llOO-l-OOll sector at Li = 8 (lower) for N = 16, Ns — 32. The origin of AKa is chosen to match the thin-torus 
ground state. The blue squares represent numerically obtained data. The assigned edge modes are labeled by dots with different colors 
corresponding to different edge environments as described in the text. The combination of two identical edge environments is marked by the 
crosses with the color of that edge environment. The combination of two different edge environments is marked by the squares filled with 
two colors, the left (right) one of which corresponds to the edge environment on the left (right) edge of the subsystem A. Here we do not 
differentiate the edge environments J^p to Ip from J-p to Ip because the former ones are just the mirror symmetries of the latter ones and have 
a momentum shift ±Ns/2. 



In the 01 IOh-1001 sector, first we start our analysis from the 
term 01100110|0110011001100110|01100110. We find two 
edge environments: 



Vp : 01100110101100110, AiV 



•Dp 



0; 



: 01001111|00100110, AiVJ = 1; 
01100100|11110010, AiV^^ = -1, 



whose counting rules are expected as 

Vp: 1,1,3,5,10, ••• ; [1,1, 3, 3, 2,1]; 
fp: 1,2,4,7,13,--- ;[1,2,1]. 



If we start from the other term, 
1001100111001100110011001110011001, we also find 
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the following two edge environments: 

V'^ : 10011001|10011001, ATVj^ = 0; 

: 10011000|11110001, ATvf = -1; 

10001111|00011001, ATvf = 1, 

with counting rules being 

V'^: 1,2,4,7,13,- ••; [1,2,4,3,1]; 
1,1,3,5,10,- ••; [1,1,2,1]. 

The combination of these edge environments can generate 
V^T)^, £p£p, V^Vp and fpfp towers in the Na = 8 sector; 
Pf'S'f, fpI^F, I^F^F ^Y^Y towers in the Na — 7 sector; 
and £p£p and £p5p towers in the Na = 6 sector. 

In the 0011+1100 sector, if we start from the term 
001 10011 [0011001100110011 1 001 10011, we find the follow- 
ing four edge environments: 

Jp : 001 1001 1 1 001 10011, A A^J'' = 0; 

: 00110010(11010011, AA^^"- ^ -1; 
-Hf : 00110000(11111100, AAr;^'^ = -2; 
Zf : 00011111(00010011, AA^J = 1, 

with counting rules being 

J-f: 1,1,3,5,10,--- ; [1,1,3,3,1]; 
^f: 1,2,4,7,13,--- ;[1, 2,3, 2,1]; 
Hf: 1,1,3,5,10,--- ;[1]; 
Xf: 1,2,4,7,13,--- ;[1]. 

The combination of edges generates towers in each Na sec- 
tor: J^F-^F and QpQp towers in the A'^ = 8 sector, J^f^f. GpHp 
and XpJ'p towers in the A^^ = 7 sector, and TpHp and If^f 
towers in the A^^ = 6 sector If we start from the other term, 
11001100(1100110011001100(11001100, we find edge envi- 
ronments Jp to Zp which are just the mirror symmetries of 
J^F to Ip with the same counting rules. For example, the edge 
environment J"^ is identified as 11001100(11001100, which 
is the mirror symmetry of J'p. Moreover, the combination 
of edge environments Jp to Ip forms towers in different A^^ 
sectors. For example, in the A^^ = 7 sector we predict and 
observe the GpJ^p tower as the mirror symmetry of the J^f^f 
tower. 



IV. ENTANGLEMENT ENTROPY 

The entanglement entropy of the Moore-Read state was 
studied earlier on the sphere and disk geometry^", and the 
topological part, 7, has been reported to be consistent, al- 
beit not in perfect agreement, with the theoretical predictions. 
However, the limitations of these geometrical setups make it 
very hard to verify whether the scaling regime (2) is reached 
or if the approximate agreement with theory is accidental. In 
addition, there are large finite size effects on the disk due to 
the (large) physical edge of the system. Here we revisit this is- 
sue using the torus setup that allows for superior control of the 



entanglement scaling properties as demonstrated for Abelian 
FQH states in Ref. 21. This method of partitioning implies 
two disjoint edges between the blocks, each of length Li, so 
the entanglement entropy should satisfy the following specific 
scaling relation 

SA-2aL,~2j + Oil/L,), 

where 7 is the topological entropy whose theoretical value is 
ln(\/4) for bosons and ln(V8) for fermions'"-". (See Ref. 51 
for the scaling relation of entanglement entropy in various 
physical systems.) 

Figs. 7 and 8 show the entropy Sa, its derivative dSA/dLi, 
and the intercept of its linear approximation, Sa — Li x 
dSA/dLi, as functions of Li in different sectors for different 
system sizes. Arguably the boson results (Fig. 7) look more 
promising. In this case, the entropy in the different sectors dif- 
fers for small Li, as can be expected from the thin torus limit 
where the 20 ± 02 states have an entropy of In 2 while the 11 
sector has zero entropy. However, from Li w 7 the entropies, 
Sa, in the three sectors are very similar (left panel), although 
the more sensitive indicators dSA/dLi (middle panel) and in 
particular Sa — Li x dSA/dLi (right panel) show some dif- 
ferences. The density entropy in the bosonic More-Read state 
appears to be about a w 0.25. In the case of fermions (Fig. 
8), we find that the scaling regime is not yet reached, even 
though one may make a crude estimate of the entropy density, 
a w 0.2. 

The entropy density of a state is an indicator of how chal- 
lenging it is to simulate the state on a classical computer, 
through a one-dimensional algorithm such as DMRG^^'^'^^, 
which has already been applied to the FQHE problem^"*, or 
through recently-proposed true two-dimensional algorithms 
like PEPS^^ or MERA^*. The larger entropy densities of 
Moore-Read states imply that they are more difficult to simu- 
late than the Laughlin states. 

Even for our largest system sizes, where we have obtained 
data for a range of ii values (A^ = 14 for bosons and A^ = 16 
for fermions), we cannot extract a reliable topological entropy 
(see Figs. 7 and 8). However, we can observe some interest- 
ing phenomena. First, the entropy densities a = dSA/d{2Li) 
of Moore-Read states are significantly larger than that of the 
fermionic Laughlin state atiy = 1 /3^' and the bosonic Laugh- 
lin state at I' = 1/2^^. Second, the entropy properties of 
bosonic Moore-Read states in the 1 1 sector and those in the 
02+20 sector become similar at large Li. Their curves of Sa, 
dSA/dLi and Sa — Li x dSA/dLi overlap after Li w 10. 

Let us also highlight some of the finite size features. For 
small Li the finite size convergence is essentially perfect and 
the curves for different system sizes are on top of each other 
in a given sector. At larger Li, the curves show a stronger 
dependence on Ng. The A^g-dependence shows up first for the 
smallest system size and at increasing Li for progressively 
larger system sizes. This reflects the fact that, for any finite- 
size system, at very large Li the edges of block A get too 
close and cannot be thought of as independent. In particular, 
once Li exceeds some value, we enter the thick- torus limit, 
and the entanglement entropy goes to some (A^ dependent) 
saturation value. Corresponding to the saturation of Sa, the 
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Li L, L, 

FIG. 7. (color online) The bosonic Moore-Read state entropy Sa, its derivative dSA/dLi and the intercept of its linear approximation 
Sa — Li X dSA/dLi as functions of Li for TV = TVs = 12 and A*' = TVs = 14 in the 1 1 sector, 02+20 sector and 02-20 sector. The theoretical 
value of the topological entropy 27 — — In 4 is indicated by the black line. In a rather large window of Li, the entropy properties of the 
Moore-Read states in the 1 1 sector and 02+20 sector are quite similar. 



derivative dSA/dLi drops to zero after some Li. Therefore, 
the appropriate scaling regime of the entropy, Sa, may be ex- 
pected to be valid only in a window of Li, after the ©(l/Li) 
term is small enough but before Sa saturates. This analy- 
sis was shown to provide excellent results for abelian FQH 
states in Ref. 21. However, as already mentioned, we find 
that the finite size corrections to the scaling are too large to 
faithfully determine the topological part, 7, of the entropy for 
the Moore-Read state. Given the limitations encountered also 
in other geometries, we conclude that an accurate and reliable 
determination 7 for the Moore-Read state remains a challenge 
for the future. 



V. DISCUSSION 

We have investigated the entanglement spectrum (ES) and 
the von Neumann entropy of bosonic and fermionic Moore- 
Read states on the torus. The ES on the torus is much more 
intricate and the analysis thereof poses a number of challenges 
compared to the sphere geometry where there is a single edge 
and a unique ground state. One such challenge is that in a 
given particle number sector, several towers appear due to pos- 
sible compensating charge transfer across the two boundaries. 

The study of the entanglement in this geometry is neverthe- 
less well motivated as it provides new insights, for instance 
by connecting to the vicinity of the microscopically well un- 
derstood thin-torus limit, and also because it may provide 
guidance for future studies of entanglement in other many- 
body systems where no natural analog to the quantum Hall 
sphere exists. In particular, the recently suggested fractional 
Chern insulators*'"'*^ are most naturally studied using periodic 
boundary conditions, i.e. on a torus. 

In this work we have suggested a procedure in order to re- 
solve the problem of the non-trivial ground-state degeneracy 
on the torus: we used exact diagonalization and choose to cal- 
culate the entanglement in the pure (simultaneous) eigenstates 
of H, Ti and r|. This is different from the mixed-state recipe 
of Refs. 33 and 42, for which we expect a superimposed en- 
tanglement spectrum and a shifted prediction of the topologi- 



cal entropy 7. 

For the ES, we found a tower structure similar to, but sig- 
nificantly richer than, what was found earlier in the ES of the 
Laughlin state. We used two complementary ideas in order 
to disentangle the ES by extending the results of Ref. 28 to 
non-Abelian states. The first approach is based on a combi- 
nation of two chiral CFT edges. Each of these is individu- 
ally similar to the edge spectrum previously extracted from 
ES studies on the sphere. This interpretation is powerful as it 
reproduces the entire ES through the assignment of a few lev- 
els. It also reflects the intricate structure of the correlations in 
the Moore-Read state: Even for one cut of our system, edges 
corresponding to different topological sectors with different 
counting rules combine to form towers. Our second approach 
uses the adiabatic connection to the thin-torus limit: A per- 
turbation analysis away from the thin-torus states yields the 
locations of the towers, and the counting rule of each edge en- 
vironment follows from a generalized exclusion principle in 
the occupation number basis. 

A further difficulty encountered when disentangeling the 
torus ES is the non-monotonic dispersion that appears for 
fermions at large Li, as the lack of a natural vacuum level 
at the bottom/center of each tower severely increases the diffi- 
culty of the assignment of the edge modes. In the present case, 
this difficulty can be circumvented by following the smooth 
dependence of the edge levels to the small Li regime, where 
the dispersion is always monotonic. 

For the von Neumann entropy, we found that the area-law 
entropy density, a = dSA/dLi w 0.20 — 0.25 (per magnetic 
length), is larger than in the Laughlin states for both bosons 
and fermions. However, the comparable smallness of a is nev- 
ertheless encouraging regarding the possibilities of simulating 
the Moore-Read state using entanglement-based algorithms. 
Our results also show that an accurate and reliable determina- 
tion 7 for the Moore-Read state on the torus remains a chal- 
lenge for the future. It is likely that alternative methods, such 
as DMRG in a cylinder setup, will be needed to reach this 
goal. 

The generalization of the analysis given here for the Moore- 
Read state to more generic non-abelian FQH states should be 
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FIG. 8. (color online) The fermionic Moore-Read state entropy Sa, its derivative dSA/dLi and the intercept of its hnear approximation 
Sa - Li X dSA/dLi as functions of Li for iV^ = 28 and Ns = 32 in the 0101 sector, 0110+1001 sector and 0110-1001 sector. The 
theoretical value of the topological entropy 27 = — ln8 is indicated by the black line. For the cut in the 0011 ± 1100 sector is equivalent 
to a translation of the cut in the 0110 ± 1001 sector. Therefore we make an average over their entropies and only show the averaged results 
(referred to as the 0110 ± 1001 sector in above). 



straight forward, but nevertheless interesting. The generaliza- 
tion to fractional Chern insulators is more challenging, but is 
likely to be rewarding. 
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Appendix A: Hamiltonian generating Moore-Read states 

The bosonic and fermionic Moore-Read states on the torus 
are the unique zero-energy ground states of translational in- 
variant three-body interaction Hamiltonians 

{k} 

(Al) 

where {fc} — fci, ^2, fcs, A;4, ^5, fcg, ak {a\) annihilates (cre- 
ates) a boson or a fermion in the state ^jjk in Eq. (2), 

+00 

E 5's,M-k/s.M~kP{{^}.{t}) 

{6},{t} = -00 

X exp I - -jT^^i + ■'^2 + ■S1S2) - -^(q + <2 + ^1*2) I 
X exp|^ti(2fc3-2fci+2si+S2)} 
X exp|^t2(2fc3-2fc2+si+2s2)|, 



and 6' is the periodic Kronecker delta function with period 
Ns- P({s}, {i}) is a certain polynomial of si, S2, ti, i2, exact 
form of which depends on the targeted filling fraction. 

We use exact diagonalization to obtain the ground states 
of (Al) after choosing a proper form of P. Up to a con- 
stant factor, when P — 1, (Al) can generate the bosonic 
Moore-Read states at filling factor v = 1, while when P = 
-i4n^n4/Ll + tl/Ll)[{s, s,f/Ll + {t, t,f / Ll], 
(Al) can generate the fermionic Moore-Read states at filling 
factor V = \/2. 



Appendix B: Edge excitation of Moore-Read state 

Compared with the Laughlin state, the edge excitations of 
the Moore-Read state are richer. It has one branch of free 
bosons and one branch of Majorana fermions obeying either 
antiperiodic (Bl) or periodic boundary conditions (B2). 

For free bosons plus antiperiodic Majorana fermions, the 
excitation spectrum is described by the Hamiltonian 

HZ^. = E [Eb{m)hlb^ + Ef{m - l/2)c|„_i/2C™-i/2], 

m>0 

(Bl) 

where b and W ic and c)) are standard boson (fermion) cre- 
ation and annihilation operators, Et,{m) [Ef{m)] is the dis- 
persion relation of bosons (fermions) and the total momen- 
tum operator is defined as K = X](Ti>o["^^»Ti^nj + (™ ^ 
l/2)c|^ -^^2Cm-i/2]- The counting rule of the edge excita- 
tions, namely the number of energy levels at each K, de- 
pends on the parity of the number of fermions (— 1)^,F = 
J2rn>o^ln-i/2^^ri-i/2- For cvcn F, the counting rule is 
1, 1, 3, 5, 10, • • • at Afc = 0, 1, 2, 3, 4, • • • ; while for odd F, 
the counting rule is 1, 2, 4, 7, 13, • • • at AA: = 0, 1, 2, 3, 4, • • • . 
Here Ak is defined as K — Kq where Kq is the lowest mo- 
mentum (Kq = for even F and Kq = 1/2 for odd F). 
For free bosons plus periodic Majorana fermions, the edge 
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TABLE I. In this table, we analyze the counting rule of the edge environment in the 1 1 sector, which is 1, 2, 4, 8, 14, ■ • ■ at Afc = 
0,1, 2, 3, 4,.-.. 



Afc = Afc = 1 


Afc = 2 


Afc = 3 


Afc = 4 




1111111111|0000 1111111110|1000 


lllllllllOlOlOO 


lllllllllOjOOlO 


1111111110 


0001 


1111111102 0000 


1111111101 1000 


1111111101 0100 


1111111101 


0010 




1111110202 0000 


1111110201 1000 


1111110201 


0100 




1111111020 1000 


1111111020 0100 


1111111020 


0010 






1111111100 2000 


1111111100 


1100 






1111020202 0000 


1111020201 


1000 






1111102020 1000 


1111102020 


0100 






1111111011 1000 


1111111011 


0100 








1111110200 


2000 








1111110111 


1000 








1111111010 


2000 








1102020202 


0000 








1111102011 


1000 








1110202020 


1000 



TABLE II. In this table, we analyze the counting rule of the edge environment Db in the 02+20 sector, which is 1, 1, 3, 5, 10, ■ at Afc = 
0,1,2,3,4,---. 



Afc = 


Afc = 1 


Afc = 2 


Afc = 3 


Afc = 4 




020202020210000 


02020202011 1000 


0202020201 10100 


020202020110010 


0202020201 


0001 






02020202001 2000 


0202020200 1100 


0202020200 


1010 






0202020111 1000 


0202020111 0100 


0202020111 


0010 








0202020110 2000 


0202020110 


1100 








0202011111 1000 


0202011111 


0100 










0202020200 


0200 










0202020102 


0100 










0202020020 


2000 










0202011110 


2000 










0201111111 


1000 



excitation Hamiltonian is 

^fdge = E [Ebim)blb,„ + Ef{m - l)4_iC„,_i], (B2) 

m>0 

for which the total momentum \s, K ~ X]m>o["^^m^m + 
(to — l)c|jj_j^c,„_i]. Through a similar analysis with that 
for the antiperiodic case, one can get that the counting rule 
is 1, 2, 4, 8, 14, • • • at Afc = 0, 1, 2, 3, 4, • • • for both even and 

oddF = Em>ocL-lCm-l- 

The counting of each edge environment observed in our ES 



should be consistent with one of the four sectors here before 
the finite size effect truncates the series after some Afc de- 
pending on the system size. 

Appendix C: The counting rules of edge environments 

Here we analyze the counting rules of the edge environ- 
ment A-R in the 1 1 sector and the edge environment I?b in the 
02H-20 sector The results are obtained by applying the gen- 
eralized exclusion rule on their thin-torus limit. All possible 
edge excitations at each Afc are listed in Tables I and II. 
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